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SUMMARY 


An instrumentation system is described which employs a minicomputer, a 
one-third octave band analyzer, and a time code/tape search unit for the 
automatic control and analysis of third octave data. With this system the 
information necessary for data adjustment is formatted in such a way as to 
eliminate much operator interface, thereby substantially reducing the 
probability for error. A description of a program for the calculation of 
effective perceived noise level (EPNL) from aircraft noise data is included 
as an example of how this system can be used. 


INTRODUCTION 

This paper describes an instrumentation system which was evolved for 
the purpose of developing more efficient and accurate dynamic data acquisition 
and analysis procedures. The instrument system utilizes a minicomputer to 
control and analyze data obtained from an analog tape recorder and a one-third 
octave analyzer. One of the major sources of error in the analysis of data 
is the improper logging of the gain values used in the recording system. A 
unique feature of this system is the utilization of automatic digital coding 
of information which, when impressed on the time code frame, is recorded 
along with the data. The recording system gain values and other pertinent 
information are then available on the tape without reference to a handwritten 
log sheet. With this information and the ability to control motions of the 
tape recorder, the analysis operation requires very few operator inputs. 

Acoustic noise data often requires additional weighting and computation 
procedures beyond the one-third octave filtering. In tie analysis of air- 
craft noise data many of these procedures are necessary to achieve a value 
known as the effective perceived noise level (EPNL). A program for computing 
EPNL has been written, providing a means of determining which factors most 
influence acoustic measurement values. This program, demonstrating all the 
features of the described computer-controlled system, will be discussed. 


DESCRIPTION OF THE INSTRUMENTATION 
General 

This instrumentation system consists of a minicomputer with peripheral 
equipment, a third-octave analyzer, a time code/tape search unit, and an 
analog tape recorder. The tape search unit, which has access to the data 
stored within the time code frame, receives command information from the 
computer, interprets the commands, and controls the tape recorder. The 


third-octave analyzer receives the analog data from the tape recorder, 
processes it, and transfers the data to the computer for further analysis. 

The heart of the system is the minicomputer to which all the other 
instruments are directly or indirectly interfaced. Some of the interfaces 
are standard; however, the interfaces to the third-octave analyzer and the 
tape search unit were not available from the manufacturers and had to be 
developed in-house speci fically for this system. 


Minicomputer 

The controlling element of the system is a minicomputer utilizing 
l6-bit words and containing 32K words of memory. The computer functions in 
a disc operating system configuration which allows for program and data 
storage on a moving head magnetic disc. Programs can be developed in 
Assembly, FORTRAN, FORTRAN IV, or Algol languages. Communication with the 
system is via a CRT terminal. Data can be output to the disc or a digital 
magnetic tape for storage, while hard copy listing and graphics are available 
on an electrostatic printer/plotter. 

The disc operating system allows segmented programs where long programs, 
divided into segments, are stored on the disc, and only one segment of the 
program is present in the computer at a time. Analyses which require many 
computational procedures or programs which provide for a number of optional 
procedures need such a capability. 1116 EPNL program which will be described 
later utilizes this technique. 

Figure 2 shows the computer and its standard peripherals (CRT terminal, 
disc, magnetic tape, and the printer/plotter ) . 


Analog Tape Recorder 

The analog tape recorder in this system can be either a T- or lU-channel, 
intermediate or wideband g: up I, FM machine (fig. 3). One channel is 
dedicated to IRIG-B time code. These tape recorders have been modified so 
that they operate in a 120 IPS, I'M mode when "tape search" is selected. 

This allows *he time code to be recorded in either FM or direct mode and 
still remain under tape search control. The tape transport of either 
machine can be remotely controlled by the computer depending on which one 
is connected to the tape search unit. Further discussion on this control 
will occur in the section on the tape search unit. Playback speeds from 
1-7/8 to 120 IPS are available. Remote control of speed selection is 
possible, but is not being employed at this time. Edge track voice annotation 
is available or an external speaker system can be utilized when voice 
annotation is on a data channel. The tape speed can be servo or tach meter 
controlled. 
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One-Third Octave Analyzer 


The analyzer used to provide the one-third octave spectral data is a 
multifilter, multiplexed detector instrument (fig. 4 ) . With this analyzer 
the spectrum can be divided into as many as 30 one-third octave bands 
ranging from 25 Hz to 20 kH... Also an integration time ranging from 1/8 sec 
up to 32 sec can be selected. The filter output values are visable via a 
nixie readout display or in bar graph form on a CRT display. The filter 
outputs are also available in parallel BCD digital form which is the output 
that is used in the interface to the computer Detailed information on this 
interface can be found in reference 3. 


Time Code Reader 

The recording of a serial time code has been in practice for some time 
and has provided a means of correlating multiple records of a single event. 
Along with the time of day, either the event number or the day of the year 
was recorded on the code. Recently the capability has been provided for 
superimposing additional information on the time code signal. As shown in 
figure 5, less than half of the serial bits in the code are needed for the 
time and event (or day) information. If the redundant recording of time 
in binary form is omitted, 1*5 bits of information can be recorded each 
second (IRIG B). 

The time code reader in this system has the control bit translation 
feature where the state of U5 control bits is displayed on LED's and is 
also provided as a parallel digital output each second. The time information 
is also available in parallel form, and both sets of digital lines are 
monitored by the tape search unit. When requested by the computer the time 
or control bit information is relayed from the tape search unit. 


Tape Search Unit 

When used in conjunction with the time code reader, a tape search 
unit (TSU) provides for considerable versatility in tape recorder operations. 
If the optional characteristic of a computer interface is provided, this 
versatility expands to an almost unlimited degree. The system described in 
this paper has a unit with these features. Figure 6 shows the time code 
reader and the tape search unit. The various operational modes will be 
described in two parts. The first, TSU-Remote, will cover operation of the 
tape recorder from the front panel of the tape search unit. The second, 
Computer-Remote, will describe the features provided by the TSU-computer 
interface. 
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TSU-Remote 


This mode is the standard TSU-tape recorder operation. All commands 
ure implemented from the front panel of the TSU. 'I*he major features of this 
mode of operation are as follows: 

(1) Set "start" and "stop" times with rotary thumbwheel switches. 

(2) "Search to Start" - The tape recorder is put into a search speed 
from the front panel "Forward" or "Reverse" pushbutton. The search speed is 
a controlled speed which is higher than the normal run speed, but it is 
lower than the wind or rewind speeds. The time code is read and the motion 
direction is corrected if improper. The tape is positioned to a point at 
least three time frames in front of the selected start time, and the 
transport is then stopped. 

(3) "Single Cycle" - The same procedure occurs as in "Search to Start" 
except that after stopping at the start time the transport is put into the 
"Drive" forward mode. An "interval" signal is provided when the start time 
is reached and is removed and the transport stopped when the stop time is 
reached. 

{h) "Recycle" - The same procedure occurs as in "Single Cycle" except 
that instead of stopping when the stop time is reached the tape is rewound 
to the start time. This cycle is continuously repeated until another 
command is given. 

(5) In addition to these automatic modes all the normal transport 
modes can be manually selected from the TSU. 


Computer-Remote 

All the features cited in the previous section are also available under 
computer control. Logic circuits in the search unit interpret commands from 
the compute. . These commands can be for tape transport control , setting 
"start" or "stop" times, or requests for time code data. The time code 
requests can be for either time or the control bit information. Documentation 
concerning this interface can be found in reference 2. 


SYSTEM APPLICATION 

With the tape recorder, the tape search unit, and the one-third octave 
analyzer under computer control, programs con be written which provide data 
analysis with a minimum of operator interface. The EPNL program description 
which follows illustrates the pertinent features of this instrumentation 

system. 


EPNL COMPUTATION - AN APPLICATION 


This section presents a description of a data analysis program utilizing 
the instrument system described in the preceding sections. The program uses 
the physical properties of measured flyover history to determine the noise 
evaluation quantity designated as effective perceived noise level, EPNL. 
Calculated in accordance with the procedures as defined in Federal Aviation 
Regulation, Part 36 (ref. l), EPNL is a single number measure of complex 
aircraft flyover noise as related to subjective response. It is calculated 
from a time sequence of tone adjusted perceived noise levels which are cal- 
culated from one-third octave bund spectra. The tone adjustments are 
determined by a procedure which estimates the extent of spectral irregularities 
at any instant of time. 


General Procedure 

As described in reference 1, the EPNL calculation is basically a 
five-point procedure: 

(1) Calculation of perceived noise levels, FNL(K), from one-third octave 
band levels determined at the K-th interval of time, where K is a number 
indicating the number of 0.5 sec intervals elapsed from an arbitrary 
reference time of the flyover signal. 

(2) Calculation of a tone correction factor, C ( K ) 

(3) Calculation of tone corrected perceived noise levels, FNLT(K), 
and determination of the maximum value, PNLTM 

(U) Calculation of a duration correction, D 

(5) Calculation of effective perceived noise level, EPNL, which is the 
algebraic sum of the maximum tone corrected perceived noise level and the 
duration correction factor. 

EPNL = PNLTM + D 


The Program 

The EPNL program was edited, loaded, compiled, and stored using the 
disc operating system of the minicomputer system. 

Many computer programs exist for the calculation of EPNL; however, this 
particular EPNL program has been written to provide various data analysis 
techniques by allowing for versatile modes of examining one set of data. 
Recause of its complexity and length, the program is structured into a main 
program and two segmented programs, operating together as a single program 
in core. 
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The main program contains all input parameters, standard and optional. 
It allows for the following operate* inputs: 

:,r_ . . 

(l) Test statistics - test name, run number, tape recorder channe 

number 


(2) Correction options - pink noise, microphone, pseudotone, ambient 




noise 


(3) Time code information - start and stop times of the calibrate 
signal, pink noise, ambient noise, and flyover data 

(U) Selection of the number of pink noise, ambient noise, and data 


spectra 







(5) Selection of the data averaging technique - linear, exponential, 







or none. 

With these inputs and the associated instrumentation, the data analysis 
program proceeds with little or no additional operator control. Similarly, 
with corresponding time code inputs, correction data can be input after 
being automatically located on the tape. Information such as amplifier gain 
settings can be read and decoded automatically from time code control bit 
information. Subsequently, the program inputs one-third octave flyover data 
each 0.5 sec increment of time from the third octave analyzer and converts 
the data from BCD (binary coded decimal) to floating point. To calculate 
an EPNL value, flyover data from 1 up to 100 spectra can be used. After 
corrections are applied to this flyover data, control is transferred to the 
first segmented program. 

The first segmented program does the actual calculation of EPNL as 
defined in the general procedure. The second segmented program controls 
the hard copy printout of the test statistics and other important parameters, 
In addition to the EPNL value, the following is output: 

(1) The spectral data printed in tabular form 

(2) A time history plot of PNLT values 

(3) Plots ol selected one-third octave spectra 

After the output is completed, control is returned to the main program 
with the operator option of additional calculations or ending of the program. 




With the segmented program feature the EPNL program is flexible, allowing 
for future amendments and additions. For example, corrections to the EPNL 
number for test day meteorological conditions and aircraft position will be 
incorporated when the instrumentation is completed for recording such data. 







A listing of the present program is contained in the appendix. 


. 




CONCLUDING REMARKS 


A computer-controlled instrumentation system to provide more efficient 
and accurate dynamic data acquisition and analysis procedures has been 
developed and an application demonstrating its use has been described. One 
of the outstanding features of this system is the ability to utilize K0 
digital information which has been coded within the time code signal. Als 
of significance is the tape search unit computer interface which provides 
for considerable versatility in remote control of tape recorder operations. 
With thi i instrumentation many of the probable physical and mental errors 
caused oy human intervention have been eliminated, thus maintaining more 
consistent accuracy in data input and analysis. In addition to minimizing 
these errors, this automation hac resulted in significant time savings for 
analysis proct uures where many channels of similar type data are being 
processed. Future developments in increased computer control could virtually 
eliminate all human error involved in processing analog data stored on 
magnetic tape 
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EPNL PROGRAM 


ORIGINAL rAU 






0301 

0092 

0003 

0034 

0005 

0036 

0007 

0G00 

0009 

0010 
0011 
0012 

0013 

0014 

0015 

0016 
0017 
0010 

0019 

0020 
0021 
0022 

0023 

0024 
C025 
0326 
0027 
0020 

0029 

0030 

0031 

0032 

0033 

0034 

0035 

0036 

0037 

0038 

0039 
CO40 

0041 

0042 

0043 
0344 

0045 

0046 
G347 
0049 

0049 

0050 

0051 


FTTI4.L 

C EPNL COMPUTATION 
C 






a -.Vnr^ „ I ,U' M C B 



c 
c 
c 
c 

C 
C 
C 
C 
C 
C 

c 

C ******.•(<* 

c 


THIS PROGRAM INPUTS 1/3RD OCTAVE DATA FROM A GR1921 ANALYZER VIA 
CALL EXEC <1.7, * ) AND CON' CRTS THE DATA FROM BCD TO FLOAT- 

ING POINT VIA CALL DCDFP. FLYOVER DATA UP TO 100 SPECTRA CAN 
USED TO CALCULATE AN EPNL VALUE. It! ADDITION TO THE EPNL VALUE* 
THE FOLLOUING IS OUTPUT* 


1) 

2) 

3) 


THE SPECTRAL DATA PRINTED IN TADULAR FORM 
A TIME HISTORY PLOT OF PNLT VALUES 
PLOTS OF SELECTED 1/3RD OCTAVE SPECTRA 



PROGRAM EPHL3.3 

DIMENSION IBFD(240 ! " * IBFA(400) # IBUF(4C0) * INAM(3) » ICD(4) 

D I f ill IS I ON SPDO C 24) . i .US ( 4) . SPP (24* 20) * SPAM(24, 20 ) . SPPC (24) 

P l MENS ION CAI.(24) ,CALA( 12) ,GD( 12) *ODA( 12) *BC( 12) . IBC<24) 
COMMON f.PL(24) * S(24) ,SN (24) * SPL 1 (24) , EPL2 (24) . SM3(24) . SMC (24) 
COMMON SPLA(24) «SPLB(24) »SPLC(24) « IC(32) « ID(62) 

common SP.V24. !2).pc(24* i2)*nc(24 * i2)*p:il(ioo)*pnlt(100) 

LG Mu i0! I on ( ICO) . SPD (24* 1U0) * OAC (24) . OAA ( 100) * NCH 
COMMON UAM. KO* PHLM. KC.OAAM. KOA* KAV 
COMMON DAV(24. 100) . ATEST(H) . I0ND(24) 

COhl 3,1 SPM(24) *KP«N« I6US(4) .KTO.KRUN 
E0U1 VALENCE ( IBFD( 1 ) -SPD( 1, l) ) 

FOU I VALENCE ( JBUF ( 1) *SPP( 1* 1)) 

LOU 1 VALENCE ( I0FA( 1) . SPAM( 1 * 1) ) 




•A* UT. CORRECTION DATA (OAC) AND 1/3RD OCTAVE BAND NOS. (IBAND) 



oncu)— 30.2 
00C(2>— 26.2 
0AC(3)— 22.5 
0AC(4) — 19. 1 
0AC(5) — 16.1 
0AC(6> — 13.4 
0AC(7) — 10.9 
OAC (B)— 8.6 
OAC (9)— 6.6 
OAC(IO)— 4.0 
OAC(ll)— 3.2 
0AC(12) — 1.9 
one (13 )— .0 
one (14) - 0 . 
0AC(l5)-.6 
OAC ( 16) ■ 1 . 
0AC(17)-1.2 







RpyROI/UCIBIUTY OF THL 
ORIGINAL PAQK IS POOR 


0052 


OOC( 18) *1 .3 

0053 


00CC19)-1.2 

0054 


OAC(20)-I. 

0055 


0f)C(21)-.5 

0056 


one (22) — . 1 

0057 


OiOC (23) — 1.1 

0058 


one (24)— 2.5 

0059 


1 Bt IP < 1 ) ■ 17 

OJG0 


1BIID(2)-I8 

0061 


1CMD(3) B 19 

C062 


1DMI)(4) "20 

0J63 


IUNDC5) u 21 

0064 


I0*ll)(6) “22 

0065 


1CHP(7)“23 

0066 


lCNP(8)-24 

0067 


ir:iD(9) -25 

0068 


irflDC 10) -26 

0069 


10IID( 1 1)-27 

0070 


IBMP( 12) -28 

0071 


irripc 13) -29 

0J72 


ICMP ( 14) -30 

0073 


IOIID ( 15) -3 1 

0074 


in:ip( 16) -32 

0075 


irriD c 17) -33 

0076 


U;;ID( 18) -34 

( 177 


IDMD( 19) -35 

0070 


10IID(20) -36 

0079 


10ND(2P-37 

COOO 


101 ID (22) -38 

0031 


10MPC23) -39 

0032 


lGMP(2.')-40 

0033 

C 


0034 

c 

HOY TOOLE CONSTOHTS 

0035 

c 


0006 


SM8(l) -0.043470 

0037 


SMO (2) -.040570 

0008 


sra (3) -.036931 

0039 


SM3( 4) -.036031 

0090 


Sill) (5) -.035336 

0091 


SrtJ (6) -.03333.1 

0092 


sra (7) -.033333 

1 1 : 


SI )(8) a .Q32051 

0034 


sra (9) -.030675 

0095 


DO 2 1-10.22 

0096 


2 srad)-o. 

0357 


sra (23) -.042285 

0098 


sra (24) -.042285 

0099 


DO 4 1-1.15 

otoo 


4 SI .C( I >-.030103 

0101 


TO 6 I t *16.24 

0102 


6 SPIC ( I 029960 


0103 

SPLAC 1)-91. 

0104 

SPLA(2) "05.9 

0105 

SPLAC 3) -07.3 

0106 

SPLAC4) "79.9 

0107 

SPLAC5W9.8 

0100 

SPLAC6) “76. 

0109 

SPLAC 7) -74. 

0110 

SPLAC9W4.9 

0111 

SPLA(9)-94.6 

01 12 

DO 8 1-10.22 

0113 

0 SPLAC I)-0. 

0114 

SPLAC23) -44.3 

0115 

SPLAC24) -50 . 7 

0116 

SPLOC 1 ) -64, 

0117 

SPLDC2) -60 . 

0119 

SPL0C3) -56. 

0119 

SPL0C4) -53. 

0120 

SPLDC5) -51 . 

0121 

SPLBC6) -48. 

0122 

SPLCC7) -46. 

0123 

SPLDCO) -44. 

0124 

SPlOC 9) -42. 

0125 

DO 9 1-10.22 

0126 

9 SPLOC I)-0. 

0127 

SPLQC23) -37. 

0120 

SPL0C24) -41 . 

0129 

SPLC Cl) -52. 

0130 

SPLCC2) -5 1 . 

0131 

SPLCC3) -49. 

0132 

SPLC C 4) -47. 

0133 

SPLC (5) •’46. 

0134 

SPLC C 6) -45. 

0135 

SPLCC7) -43. 

U 136 

SPLC CO) -42. 

0137 

SPLC ( 3) -41 . 

0139 

DO 10 1-10.14 

0139 

10 SPLCC I) -40. 

0140 

SPLC C 15) -30. 

0141 

SPLCC 16) -34. 

0142 

SPLC C 17) -32. 

0143 

SPLCC 10) -30. 

0144 

SPLCC 19) -29. 

0145 

SPLCC20) -29. 

0146 

SPLC (2!) -30. 

0147 

SPLC (22) -31 . 

0143 

SPLCC 23) -34. 

0149 

SPLC (24) -37. 

0150 C 


0151 

IHAMC 1 ) -2HCA 

0l52 

IMAMC2) -2HLC 

0153 

1MAMC3) -2H 


10 


RWKUDUCrbD.ITY Or THE 
nRK',lNA4. PAGE IS IVOR 






URITE(Ml) 

11 FORMATC "TYPE : ■/ 

C“ 1 -START OF LALCALATIONS"/ 

C" 2-ADD ITIONAL CfiLCALOTIONS*/ 

C* 3-END Or CALCALATI0N5“) 

REAPC 1 #*)MCAL 
GO T0( 14.07. 12JMCAI. 

12 STOP 

14 LIRITEU. 16) 

1A FORI ‘.AT ( “ENTER 5 OR S LETTER TEST NAME. EXAMPLE: YUMA II* > 

FEAPC 1 . 15)ATEST 

15 FORMAT (3A2) 

IIR1TEC1. 13) 

13 FORMATC “ENTER RUN NUMBER") 

READ<1.*)KRUN 

MMTCC1. 1 7) 

17 FORMATC “ENTER NUM.1ER OF DATA CHANNELS*) 

FEfiDC l.*)JIICH 
UR I TEC 1. 10) 

10 FORMATC “AMBIENT NOISE CORRECTION ? TYPE 1 “YES. 2-N0“) 

READ( l.*)KAM 
ID tfv I Tfs < 1.20) 

2u rnRMAT(“Pir::i noise correction ? type i-yes. 2»no“> 

rEnnci.aoKORP 
UR I TEC 1.22) 

22 FORI v iT < “PSEUDOTONE CORRECTION? TYPE 1- YLS. 2- N0“) 

READC 1,*)KT0 
WRITE C 1.25) 

25 t o. ;: iTCI 1 IC FPEO RESP - U1ND SCREEN CORR ? TYPE 1-YES. 2-N0“) 

RFODC 1 «*)KORM 
IF ( 1 -KORM) 35.20 

ENTER MIC. FREQ. RESPONSE AND HIND SCREEN CORRECTIONS VIA 
PAPER TAPE 


20 NCH-O. 

2D NCIMICH-H 

30 (JRITECI.31) 

31 FCMATC INSERT MIC. CORR. TAPE. TYPE :GQ") 

PAUSE 

REAP(5,*)RC 

IFCNCH-JNCII)29.36 

if r:n microphone correction, set correction parameter to zero 

35 DO 35 MCH-1, JNCH 
DO EG I -I. 24 
RC(I.NCH)-0. 

33 CONTINUE 
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0205 

0286 

0207 

0208 

0209 

0210 
0211 
0212 

0213 

0214 

0215 

0216 

0217 

0218 

0219 

0220 
0221 
0222 
0223 
P224 

0225 

0226 
0227 
0 223 

0229 

0230 

023 1 
P232 

0233 

0234 

0235 
0235 
0237 
0233 
0239 
02.10 
0?1l 

0242 

0243 

0244 

0245 
0245 
0247 
0249 

0249 

0250 

0251 

0252 

0253 

0254 

0255 


SET "JP FOR CAL. SIGNAL 

39 IJRITE(1«40) 

40 F01!MAT( "READY FOR CAL SIGNAL - ? ENTER DB VALUE FOR CAL LEVEL') 
READ< 1,*)CC 

41 WRITE ( 1.42) 

42 FORMAT (/'ENTER TIME CODE AT START (DAY. HR.. MIN., SEC.)'/) 
REAP( 1,*0 I DUS (4) , I BUS (3) . I BUS (2) . IBUS ( 1 ) 

WRITE (1.44) 

44 FORMAT(/'HNTER STOP TIME (DAY. HR.. MIN.. SEC.)'/) 

READ ( 1 , *) M3US ( 4) . MSUS (3) . M3US (2) . M8US ( 1 ) 

45 WRITE (1.46) 

46 FORMAT (/'ENTER CHANNEL NUMBER (1-12) '/) 

READ( l.*)NCH 

CALL EXECC2. 13143. IBUS. 430) 

CALL EXEC(2. 1414Q.MQUS.400) 

CALL EXEC(2. 15140.129.400) 

43 CALL EXEC(1. 1149. ICD, 1) 
i;:CH=IAND(ICD(2).740B) 

INCH" INCH/32 

IFdMCII.riE.MCH) GO 10 43 
ITG"IAND( !CD(2) , 17B) 

IUG n I AHD( ICD (2) , 200) 

IF ( IUC.NE.0) !Ul>5. 

DCdlCH) "FLOAT ( IO kITG+SO+IUG) 

CAlL E4EC(2. 15143.16.100) 

CALL EXEC (2. 15140. 130. 103) 

calacncip-o 
rn go hcm.-i.io 

CALL EXEC (1.7. ICC, 24) 
ent.L GCDFP(IDC.CAL.24. IERR) 
ir-(O-IERR) 50.54 
50 WRITE (1.52) 

52 FOnMAT(/// rt DATA OVERLOAD!") 

CALL FXEC(2. 1514B. 16, 100) 

CO TO 3' 

54 CAI.M-0 

DO 59 1“ 1 .24 
IF (CAL ( I) .GE. CALM) 56, 59 
56 CALM-CALd) 

59 CONTINUE 

CAL A ( MCI I) -( (CALM-COLA (NCH) ) /NCAL ) +CAL.A (NCH) 

60 CONTINUE 

IF (HCH- JI'JCH) 62. 64 
62 WRITE (1,63) 

G3 FORMA T » / “ RE A D Y FOR NEXT CAL SIGNAL.'/ 

C'NEIJ TIM* CODE? TYPE: 1«YE3, 2-NQV) 

READ( l.HOKTC 
IF(KTC.EQ.2)G0 TO 45 
GO TO 41 

64 IF(1-I(0RP)80,65 




, c©^\s S 00 ^ 


0256 C 

0257 C SET UP FOR PINK NOISE 

0258 C 

0259 65 WRITE (1.66) 

0260 66 FORMAT ( "SELECT NO. OF PINK NOISE SPECTRA. MAX-20') 

0261 READ ( 1 . *)NP 

0262 IF (20-IJP) 65 *68 

0263 C 

0264 C ENTER PINK NOISE 

0265 C 

025S 63 WRITE Cl >70) 

0267 70 FORMAT (/'"READY FOR PINK NOISE: ENTER TIME CODE AT START V 

0260 C" (DAY. HR.. MIN. .SEC. ) 

0269 READ ( 1 . *■) I0US<4) . ICUS<3) . IBUSC2) . IBUS( 1) 

0270 WRITE (1.72) 

0271 72 FORMAT ( /"ENTER STOP TIME (DAY. HR.. MIN.. SEC.)'x) 

0272 PEADC l.«.)MJUS(4) ,M13US(3) .M3US(2) .MSUS(l) 

0273 73 WRITE (1.74) 

0274 74 FORMAT (/"ENTER CHANNEL NUMBER (1-12)V) 

0275 PEAD( l.tfONCH 

0276 CALL EN“CC2. 13143. IBUS.4O0) 

0277 CALL EXEC (2. 14143, MBUS. 400) 

0278 CALL EXEC (.2. 15140. 130.400) 

0279 NPP»NP:is24 

O^Sa CALL F;T C( 1.7, IBUF.NPP) 

0231 CALL BC1‘ R( ICUF.SPP, MRP. IERR) 

0232 IF (0- IERR) 73. 00 

0233 73 WRITE (1.79) 

0234 79 FORMAT ( •DATA OVERLOAD!") 

0235 CALL EMEC(2. 1514B. 16. 100) 

0236 GO TO 63 
0337 C 

02G3 C ALL INITIAL VALUES OF AVERAGE PINK NOISE SPECTRUM SET TO ZERO 
0239 C 

0290 80 DO 84 1-1.24 

0291 SPA ( I . NCH) °0 . 

0292 IF ( 1-KORP) 84. 02 

0293 C 

U294 C IF PINK NOISE ENTERED. RESET PINK SPECTRUM TO AVG. OF HP SPECTRA 
0_95 C ENTERED 

0296 C 

0297 02 DO 04 K-l.NP 

0299 SP.AU.Hni) -SPA( I.NCH)+SPP( I .K) /FLOAT(NP) 

0299 84 CONTINUE 

0300 C 

0301 C CALCULATION OF PINK NOISE CORRECTION VALUES 

0302 C 

0303 DO G5 1=1.24 

0304 86 PC ( I. NCI !) -SPA (0. NCH) -SPA (I. NCH) 

0305 IF (K0RP.E0.2) GO TO 07 

0306 IF (NCH- JNCH) 33, 87 
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0307 03 WRITE (1.85) 

0338 05 FORMATC/ "READY FOR NEXT PINK NOISE CHANNEL."/ 

0309 C'NEW TIME CODE? TYPE: I -YES. 2-NO"/) 

0310 READC 1.40KTC 

0311 IF (KTC.EQ.2)G0 TO 73 

0312 GO TO 68 

0313 87 IR1-KAM) 102,80 

0314 C 

0315 C SET UP FOR AMO I ENT DATA 

0316 C 

0317 83 URITEC1.89) 

0318 09 FORNATt "SELECT NO. OF AMBIENT DATA SPECTRA CNA) . MAX-20") 

0319 READC l.*)NA 

0320 IF (20-NA) 83.90 

0321 90 UR I TEC 1.92) 

0322 92 FOHMATC "READY FOR AMBIENT DATA: ENTER TIME CODE AT ST^RT"/ 

0323 C" (DA?. HR.. Mill.. SEC.)"/) 

0324 READC 1.*) I BUS (4) . I0USC3) . I BUS (2) . I BUS ( 1 ) 

0325 WRITE ( 1.94) 

0326 94 FR^MATC/ "ENTER STOP TIME (DAY. HR.. MIN.. SEC.)"/) 

0327 F (l.*)M3US(4).NSU J l3).N i 2USC2).M3USCt> 

0320 W: '£(1.96) 

0329 96 F0RMAT(/"EHTER CHANNEL NUMBER (1-12)“/) 

0330 READC l.*)NCH 

0331 CALL EXECC2. 13143. IBUS.400' 

0332 CALL EXECC2. 14140. MBUS. 400) 

0333 CALL EXEC (2. 15 143. 130.400) 

0334 93 CALL EXEC ( 1 . 1 1413. 1CD. 1 ) 

0335 INCH- IANDCICDC2). 7400) 

0336 INCH- IHCH/32 

0337 IF C INCH. NE. MCI I ) GO TO 98 

0333 ITG- I AND ( ICD C2) . 17B) 

0339 IUG-IANDC ICDC2) .200) 

0343 IFCIUG.NE.O) IUG-5. 

0341 BDACNCH) -FLOATC 10* ITG+60HUG) 

0342 CALL EXECC2. 15149, 16. 100) 

0343 CALL E)EC(2. 15 149. 130, 100) 

0344 NAA«NA*24 

0345 CALL EXECC1.7. IBFA.HAA) 

0346 CALL BCDFPC IO. A.SPAM.NAA. IERR) 

0347 IF (0- IERR) 100.102 

0343 100 URITEC1.101) 

9349 101 FORMAT ( "DATA OVERLOAD I") 

0350 CALL EXECC2. 15140,16, 100) 

0351 GO TO 90 

0352 C 

0353 C SET AMBIENT SPECTRUM TO AVG. OF NA SPECTRA ENTERED 

0354 C 

0355 
0556 
0357 


lU 


102 DO 104 1-1.24 
SPMCD-0. 

IF ( 1-KAM) 104. 103 


KKf'Ki »hl ('{)<!/ n\ 


0358 

0559 

n 7 bO 

0361 

0362 
03t>3 
0 T b4 
0 ‘55 
0366 
0567 
C3G8 
0369 
oiro 

037 1 

0173 
(1 ' < 

0374 

0375 
0375 
LI ;/7 
0<70 

V '70 
non 

I 

0332 

0303 

IRM 

i' 

I' 35 

‘ 

0 J0 

I 

Ov'O 
0301 
0 S92 
0503 

I 

0305 

0>06 

0107 

o ;oi3 
0399 
0 "1-1 
0 ’ )| 
0132 

0 1r J3 

1 ) » 
o :o3 

0405 

V 1.17 
0403 


103 PO 101 K-l.NA 

5PM( I)“SPM(l) ♦SPAM(I.K)/FLOnT(NA) 

104 CONTINUE 

IF ( 1-KAM) 122. 105 
C 

r CORRECTION OF AMBIENT DATA FOR PINK NOISE. MIC.. CAL. LEVEL. 
C RAINS AMD ANALYZER SETTINGS 

C 

105 PO 106 1-1.24 

106 SINK I ) -SPM( I ) +PC ( I . NCH' +RC ( I . NCH) + (CC-CALA(NCH) ) 

C K 8 DA (Nf H)-BC(NCH) ) 

C 

C ENTER RUN PARArETER PATA 

C 

122 I.C?1TE( l. 123) 

123 rU *t:iT( "SELECT NO. OF PATA SPECTRA (N) . MAX-100') 

REAP< 1 .iON 

1FC100-N) 122.124 

124 iJRITECI. 125) 

125 FPMMATC 'SELECT AVG TEC .;f!UE. TAPE 1-LIN.. 2-EXP. . 3-NONE') 
REAO( 1 . *)KAV 

l?7 UMITEd. 120) 

139 FU ’MAT ( 'READY TO COMPUTE EPNL. ENTER TIME COPE AT START'/ 

C* (PAY. MR.. MIN.. SEC.)"/) 

P iH 1.*) IGUSC4) . IUUSC3) « 1BUSC2) . IEUS(l) 
u;iT2 ( l. 130) 

130 fl 3M.1T(/"EN TER STOP TIME (DAY. HR.. MIN.. SEC.)'/) 
r AD ( 1 . *) M3US ( 4) . MOUS ( 3 ) . T BUS (2) . MOUS ( 1 ) 

1.'?1TE( 1. 132) 

132 El 'I ; )T(/"LNTER CHANNEL NUMSER (1-12)"/) 

Rl- ADU.*)NCH 

CA> L E>'“C(2. 13 1.10. IOUS.40O) 

CALL h*yec;2. 1414B. MOUS. 400) 

FALL F7\“C(2. 15140. 130.400) 

134 CALL E> 3C( 1. 1 140. 1CP. 1) 

INCH- IAMD(ICD(2). 7400) 

INCH- i NCH '32 

Ir( INCH. Nil. NCH)GO TO 134 
1 1 G n 1 Pi ID ( 1CD(2) .170) 

1UG • IrtMO ( ICD(2) . 200) 

IT ( IUG.NE.0) IUG-5. 

F0(N :il>‘; LOOT l 10 ■ ITG+hd HUG) 
f HI F' C(2. 15140. 16. 100) 
l ” L E: :C(2. 15 141), 130. ICO) 

ND-N 74 

cc.l p :ct 1.7. irrn.NPi 
i ah nc?ri»( ir. p.srD.MP. ierr) 

If (U- I ERR) 155. 150 

136 U3ITC( 1. 137) 

137 n -mrcPATA ovr::LuAPi*) 

CALL EXEC (2. 15 140. 16. 100) 
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0.109 GO TO 127 

0410 C 

O-lll C CORRECTION OF RAU DATA FOR PINK NOISE. MICROPHONE. CAL. LEVEL 

0412 C GAINS AND ANALYZER SETTINGS 

0413 C 

0414 150 DO 151 K-l.N 

0415 DO 151 1-1.24 

0416 151 SPD( I .K) -SPD( I.K)+PC( I .NCH)+RC( I.NCH) + (CC-CALA(NCH) ) HBD(NCH)- 

0417 CPC(NCH)) 

0410 IF(l-KAM) 190. 150 

0419 C 

0420 C CORRECTION OF RAU DATA FOR AM0IENT NOISE 

0 421 C 

0122 159 DO 100 K-l.N 

0-123 SPDC( 1) -0. 

0424 DO 176 1-1.24 

0425 SPDC( 1+1) -0. 

0125 IF ( CSPDC I,K)-SPNfI)> .GT. 10.) 164. Iu0 

0427 160 IF((SPD( I.K)-SPM( I)) .GT.5) 162. 166 

0420 162 SPD(I.K)-1O+(ALOGC10.**.;SPD<I.K)/10.) 

0129 C-10.**(SPri(I)/lO.))/2 ■’ 

0430 164 SPDC( I)-SPD(I.K) 

0431 GO TO 176 

0432 166 IF(12-I) 172.160 

0433 169 IF(SPD(I+1.K) . EQ.SPDCC 1+1) ) 170. 176 

0434 170 SPi)a.K)-SPD(I+l.K)-(((SPD(I+2.K)-SPD(I+l.K)) 

0435 C+(SFD( 1+3.K)-SPDC I+2.K) )+(SPD( I+4.K)-SPD< 1+3. K)) 

0436 C + (SPDC I+I'.KI-SPDC I+4.K) ))/4. ) 

0437 GO TO \ 64 

C 139 172 IFiSPD(I-l.K) .EO.SPDC( I-D) 174. 176 

0439 174 S^D ( I ,«) -5PD( I-t.K)-(( (SPD ( 1-2. K) -SPD (I-l.K)) 

0440 C+(5PD( 1-3. K) -SPD ( I-2.K) )+(SPD( I-4.K)-SPD( 1-3. K)) 

0-141 C + (SPD ( 1-5. K) -SPD < I-4.K) ) )/4. ) 

0442 GO TO 164 

0443 176 CONTINUE 

0444 DO 100 1-12. 1 - - 1 

0145 IFCSPD(I.K) .EO.SPDCC I) ) ISO. 170 

U44S 170 SPD( I .K) *SPD( I+l,K)-( ( (SPD< I+2.K)-SPD< 1+1 ,K) ) 

0447 C-* (SPD( I+3.K)-GPD( 1+2.10 )+(SPDC I+4.K)-SPD( I+3.K) ) 

0443 C+(SPD( I+5.K)-SPD(I+4.K)))/4.) 

0^19 SPDC( I)-SPD(I.K) 

0:30 100 CONTINUE 

3451 190 L r ?ITEf6. lOl)OC(NCH) .BDCNCH) 

C‘32 101 F( RMAT(2X, lGliCAL GAIN- . I4/2X. 1 1HDATA GAIN- .14//) 

04J3 CULL EXEC(O.INAM) 

0454 END 
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0001 

0032 

0003 

0004 

0005 
0036 
0007 
0039 

0039 
0010 
0011 
0012 

0013 

0014 
0315 
0016 
0017 
0313 

0019 

0020 
0021 
0022 
0323 
0024 

0325 

0326 

0327 
0329 
C329 
0030 

0331 

0332 

0333 
0034 
6035 

0336 

0337 
0333 
P339 

0040 
0341 
6342 
0043 
6344 

0045 

0046 

0047 
0043 
0049 


C 

c 

c 


c 


c 

c 


c 

c 

c 

c 


c 

c 

c 

c 
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PROGRAM CALC. 5 
DIMENSION INM<3> 

COMMON SPL (24) . S (24) . SN (24) . SPL 1(24). SPL2 (24) . SMB (24) . SMC (24) 
COMMON SPLA(24).SPLD(24).SPLC(24). IC(32). ID(G2) 

COMMON SPA(24. J2).PC(24. 12).RC(24. 12) .PNL( 1C0) .PNLT( 100) 
COMMON OA ( 1G0) »SPD(24, 100) .0AC(24) ,OAA( 100) .HCH 

com ;oi i ohm. uo . pnlm. kc , oaom. koa. kav 

COMMON DAV(24, 100).DTEST(3). IBND(24) 

COMMON SPM(24).KP,N. , BUS(4) .KTO.KRUN 

INM( 1) °2HLI 
INI1(2) "2HST 
11!M(3) **2H 

PRESET EPML PARAMETERS 
OAM=0 . 

PNLM-0. 

CMAX-0. 

BEGIN CALCULATION OF 3-SPECTRUM AVG.. OVERALL. "A* UT. . 

NOY. TONE CORRECTION AMD SAVE AVG. SPECTRUM DlW(I.K) 


DO 400 K-l.N 
IF ( 1 — K ) 194. 192 

192 DO 193 1-1.24 

193 SPLCD-SPOU.K) 

GO TO 214 

194 IF(2- k'AV) 192.210. 196 
3-SPECTRUM LINEAR AVERAGING 


196 IF (2-K)204. 199 
120 DO 202 1-1.24 

230 SPL ( I ) =-3.0103+40 . +10 . *(ALOG ( 10. <•*( (SPD ( 1 . 1 ) -40 . 1/MO . ) 
C +10. *.-:<((SPD( I.2)-40.)xl0.))/2. 30259) 


202 CONTINUE 
GO TO 214 
204 DO 203 1-1.24 


2CS SPL ( I ) —4 . 77 1 2+40 . * 1 0 . * ( ALOG C 1 3 . ** ( ( SPD ( I . K-2 ) -43 . ) / 1 0 . ) 

C +10. M-:t( (SPD( I.K-D-40. )/10. ) + 10.siM:( (SPD( I.IO -43. )/10. )) 42.30259) 
209 CONTINUE 
GO TO 214 
C 

C EXPONENTIAL AVERAGING 
C 
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B050 

0051 

0052 

0053 

0054 

0055 

0056 

0057 
0053 
0053 
0060 
0061 
COS 2 
0063 
C3S4 
0065 
0366 
0067 
0069 

0069 

0070 

0071 

0072 
6973 

0074 

0075 
C 076 
C 177 

0073 
C379 
COCO 

0031 

0032 

0033 

0034 

0035 

0036 
6087 
0099 
0039 

0090 

0091 
P q 92 

0093 

0094 
CJ35 
0035 
6097 
0093 

0099 

0100 


210 DO 214 1-1.24. 

212 SFLCn»10.*ALOG(10.**(SPL<I)/10.) + (l.-EXP(-.55))* 
C ( 10. x< k(3PD( I .K)/10. )-10.**(SPL( 1)^10.) ) )/2. 30259 
214 CONTINUE 
C 

C OVERALL CALCULATION 
C 

0A(K)-5PL(1) 

DO 220 I a 2.24 

DIF = ABS(SPLU) -OA(K)> 

ADD=lO.*ni.OG(l. + lO.**(-DIF/ia.)>/2.30259 
1F(SPL( I ) -OA(IO )2 16.2 18 
216 0A(1O = COCK) +ADD 
GO TO 220 

210 OA(K) = SPL< I)+ADD 
220 CONTINUE 

IF (OAM-OA(K) )222.224 

222 OAM-OACK) 

K0-K 

C 

C "A" UT. CALCULATION 
C 

224 OAA v K ) =SPL ( 1 ) +OAC < 1 ) 

DO 230 1-2.24 

D I F -ASS ( SPL ( I) +OAC ( I ) -OAA <K) ) 

ADD- 10. a-fll.OG (1. + 10. »*C-D IF/10. ))/2. 30259 
1 F (SFL ( I > ■} CAC ( I ) -OAA CIO ) 226. 228 
226 OAACIO-OAOCO+ADD 
GO TQ 230 

223 OAACIO-SPLCn+OACm+ADD 
230 CONTINUE 

IF (OAAIi-OAA (K) ) 232. 234 

232 OAAM-UAA (K) 

KOA-K 

C 

C PNLTCIO CALCULATION 
C 

234 SUMN-0. 

SNM-0. 

1=0 

236 I-I+l 

IF (24- I ) 253.239 

233 IF (9- 1)240.246 
240 IF C22-D24S.242 

242 IF(SPL( *)-SPLC( I) ) 252.244 

244 SIIU) «E/PC2.30259’!:SMC( I)>K(SPL( I)-SPLC( 1 ) ) ) 

GO TO 254 

246 1F(SPLC1)-SPLB( I) >252,243 

243 IF (SPL ( 1 )-SPLA ( I ) ) 250, 244 

250 SM (15 -DP (2 . 30259*31*® ( I ) *<SPL C I ) -SPLD ( I) ) ) 

GO TO 254 
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0101 252 SN(I)-0. 

0102 254 SUMN-SUTIM+SNCI) 

0103 IF <SNN-SN< 1 ) >256.236 

0104 256 SNN-SN(I) 

0105 GO TO 236 

0106 253 TSN»0.G5:SNN+0. 15#SUriN 

0102 PNL(K) -4 J. 0+33. 3*(flL0G(TSN)/2. 30259) 

0109 IF(PNLM-PNL(K) )260.262 

0109 260 PNLN-PNl(K) 

0110 KP-K 

0111 C 

0112 C PNLT(K) CALCULATION (10 STEP TOME CORRECTION) 

0113 C 

PI 14 262 l>0 264 1-1.3 

0115 S(I)-0. 

0116 264 SPLKI)-SPL(I) 

0117 1-3 

0118 265 1-I+i 

0119 IF(24- 1)305.266 

0120 266 SC I ' - SPL ( 1 ) -SPL ( I - 1 ) 

0121 IF (5. -PBS <S( n-SU-D) >268.205 

0122 263 IT'SC D ) 270. 270.274 

0173 270 ircSC !-l))205.235.272 

0124 272 SPLl<l-l)>(SPL(!)*5T l l C!-2)>/2. 

0125 GO TO 285 

0126 274 !F(5( I — 1 > — S ( I ) ) 276. 205 

0127 276 IF (23- 1)273,280 

0129 27J SPLl(l)«SPL(t-l)+S(I-l) 

0129 GO TO 265 

0130 2CD SPL 1 ( I) «(SPL( 1+1 ) +SPL( 1-1 ) )/2. 

0131 GO TO 265 

0132 205 SPLKI)-SPL(I) 

0133 GO TO 265 

01<4 305 SPL2 (3) -SFL (3) 

0135 SPl 2 ( 4) -SP L2 ( 3 ) + ( SPL I ( 4) +SPL 1 ( 5 ) -2 . +SPL 1 ( 3 ) ) /3 . 

0136 1-4 

0137 310 l-l + l 

0133 IF (23- 1)320.315 

0139 315 5PL2( I ) -SPL2( I- 1 ) +(SFL I ( 1+1 ) -SPL 1 ( 1-2) )^3. 

0140 GO TO 310 

0111 37 1 SPL 2 (24) -SPL2 (23) + (.?.•: SPL 1 (24) -SPL 1 (23) -SPL1 (22) )/3 

0142 CN r 0. 

0143 1-2 

0144 325 1 *1+1 

0145 IF (24- I) 330. 330 

0146 330 F J SPL ( I )-SPL2( I) 

0147 IF (F-3 . ) 325,335 

0149 335 IF (F-2Q. )340.345 

0149 340 C-F/6. 

0150 GO TO 350 

0151 345 C-10./3. 



6152 


350 

IF (21-1)370.355 

01153 


355 

IF C 10-1)360.370 

0154 


360 

02. *C 

0155 


370 

IF (CM-C)371.325 

0156 


371 

IFU-KT0537S.372 

0157 


372 

lFi 13-1)375.325 . 

0153 


375 

Cll-C 

0159 



GO TO 325 

0160 


330 

FNLT<K)-PNL(K)+Cri 

0 161 



IF (Cnox-cri) 390.394 

016? 


390 

CMAX-CM 

0163 



KC-K 

0164 

C 



0165 

c 


SAVE AVG. SPECTRUM 

0166 

c 



0 167 


394 

PO 395 1-1.24 

0163 


395 

PAVCl.K)-SPL(I) 

0 169 


400 

CONTINUE 

0170 



CALL D'ECC8. INM) 

0171 



CALL EPNL3 

0172 



END 


■* * MO ERRORS* 
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0001 c 

0002 C 

0003 PROGRAM LIST. 5 

0004 DIMENSION SPDO (24) . SPDP (24) 

0305 COMMON SPL ( 24) . S (24) . SN (24) , SPL 1 (24) . SPL2 (24) . SMS (24) , SMC (24) 

0(36 COMMON SPLA(24) .SPL8(24) .SPLC(24) . 1C(32) . ID (62) 

Pf. 17 COMMON SPA(24. 12).PC(24. 12).PC(24. 12) .PNL( 103) .PNLT( 100) 

0339 COMMON 03 ( 101) . SPD (24. 103) . OAC (24) . 0AA ( 100) . NCH 

0039 COMMON OAM. KO. PNLH* KC. OAAM. KOA. KAV 

0010 COMMON DflV(24. 103) . ATEST(3) . IDND(24) 

OC 1 1 COMMON SPM(24) ,KP.N. IBUS(4) .KTO.KRUN 

0012 C 

0013 PNLTM-0. 

0014 KCOP-O 
OC 15 C 

0016 C LISTING OF TFST NAME. RUN NUMBER. AMD INPUT PARAMETER DATA 

0017 C 

PO 18 LI? I TE ( 6 . 435 ) ATEST. KRUN . NCH 

0019 405 FORMAT <2X. 11 HTEST NAME: . 3A2/2X. 7HRU1I NO. . I2/2X. "CH NO. •.13//) 

0C20 DO 431 I = 1.24 

CC 21 SPDO ( I ) °CPA ( I . NCH) 

0022 401 CONTINUE 

0323 DO 402 1=1.24 

0C24 SPDP ( I ) "PC ( I . NCH) 

0C25 402 CONTINUE 

0026 WRITE (G.40S)SPD0.SPDP 

0C27 405 FORMAT (91 1 PINK SP ,24( 1X.F4.2)/4X.5I!C0RR .24( 1X.F4.2)/) 

GC20 DO 407 1=1.24 

0529 SPDO ( I ) =RC ( I » NCH) 

0333 437 CONTINUE 

0031 IJRITE (6. 433) SPDO 

09.12 403 F0RMAT(9!I MIC CORR. 24( IX. F4. 2) ///) 

UG33 WRITE (6.409)SPM 

0334 409 FORMAT (9H AM3. SP .24( 1X.F4. 1)/) 

0315 C 

0036 C SPECTRA DATA TADLE 

003? C 

0333 410 IJRITE (6. 420) 

0339 420 FORMA T(2X. II IK. G! I 0ASPL.5H 50 .5N 63 .5H C3 .5H 100 .5H 125 

0040 C5H ISO . 5M 209 .511 250 .5H 315 .5H 430 .511 5C0 .5N 630 .5H 000 

CPU C5H IK .5H 1 .25.511 1.6 .5M 2.0 ,5N 2.5 .511 3.15.511 4.0 ,5H 5.0 

C942 C5! I 6.3 .5H 0.0 .5H 10K /) 

G043 DO 423 K=1.N 

CO M DO 421 1 = 1.24 

03-15 421 SPDO< I) «DAV( I.K) 

0346 425 l'2ITE(6.427)K.0A(K) .SPDO 

0947 42? FORMAT ( I4.25( 1X.F4. 1 ) ) 

0943 423 CONTINUE 

0043 C 
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0050 

0051 

0052 

0053 

0054 

0055 
0055 
0057 

0050 

0059 

0060 

0051 
0CG2 
0033 

0054 

0055 
00 55 
C 957 
0053 

0069 

0070 

0071 
C372 

0073 

0074 

0075 

0076 

0077 
0079 

00 79 
0039 

1 

CC02 
0UG3 
0° ’ 1 
I 

0056 
0097 
0 1 
CO . 3 
0090 
OC '1 
co ;2 

0033 

I 

0095 
003G 
C097 
00 , ! 9 
0099 
0103 


C DURATION CALCULATION 
C 

LT1-0 

LT11-0 

LKT-0 

TOT-O. 

DO <149 K-l.N 
IF (PliLTN-PNLT (IO )43O.440 
430 PNLTn-PNLTCK) 

KT-K 

440 CONTINUE 

DO 530 K-l.N 

lF(Pi;LT(IO + 10.-PNLTM)495.460 
450 IF ( 1-KM70. 465 
455 LT1-1 
LT1 1 - 1 
CO TO 5G0 

470 IF (LTD 433. 475. 485 
475 l.Tl-K 

IFtLTl 1)580. 490. 5S0 
430 LT11-K 
GO TO 5C0 
405 IF (K-N)5SO»490 
493 LT2-N 

GO TO 594 

495 1FCLTDL01.58G.S01 
501 LT2-K-1 

504 DO 510 M-LT1.LT2 

310 TOT-TOT+10.**( (Fl)LT(rt)-40.)/'10.) 
LK-LT2-LT1+1 
LKT "LKTh-LK 
LT1-0 

580 CONTINUE 

LKT- I AND (LKT. 1) 

IF(0-LK7)535.G05 

533 I."(PNLT(LT1 1) -PNLT (I.T2) ) 595. 590 
590 KX-LT11 

LT11-LT11+1 
GO TO 596 
595 KX-LT2 
LT2-LT2-1 

595 TOT TOT- 10.*"- ( (PNLT(KX) -43. ) 4 10. ) 

505 L'UR-in. • , '.LOG ( TOT) /2 . 30259-PNLTi 1- 1 3 . +43 . 

C 

C EPNL CALCULATION 
C 

EPNL-PNLTM+DUR 

C 

C 

C TINE HISTORY PLOT OF PNLT. EACH X- 0.5 D0 
C 
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ft 


0 10 : 
0102 

0103 

0104 

0105 
010(5 
0107 
0ing 

0109 
0113 
0111 
0112 

0113 

01 14 

0115 

0116 

01 17 

0110 

0119 

0120 
0121 
0122 

0123 

0124 

0125 

0126 
0127 
Oil 8 

0129 

0130 

0131 

0 132 
013J 

0134 

0135 

0136 

0137 

0133 
P139 

0140 

0141 

0142 

0143 

0144 

0145 

0146 

0147 
0149 

0149 

0150 

0151 


RFPKODUCDIILITY of thf 
ORIGINAL PAGE IS POOB 


ITE (6.593) ICUS(4), 1GUS(3). IBUS(2). IBUS(l) 

593 F0RNAT(2/3X. 19HSTART TIME AT Kol- .3(12. Ills). 12) 

G90 UR ITE (6.900) 

900 FOP! 1T(2/3X. 1 Hl< . 36X. 4HPNLT . 5X. 3HPNL . 6X. 5H0ASPL. 5X. 3HSLA/) 
DO 930 K-l.N 

I NX" IF IX(2.*(PNLTM-PNLT(K) ) ) 

I F ( 29- I NX) 9 1 0 . 920 
910 IMX-29 

920 DO 925 11= 1.30- 1 NX 
925 IC (ID =33 

DO 923 N=31- INX.32 
923 1C (ID =32 

UNITE (6. 930) K. IC.PIILT(K) .PNL(K) .0A(K) .OAA(K) 

930 FORMAT ( [4. 2X. 32R 1 . 4C2X, F7 . 2) ) 

C 

C IMPORTANT PARAMETERS LISTING 
C 


UR I TE ( 6 . 943 ) OAI1. ICO. PNLM. KH . CfiAX. KC . PNLTM. I'T. OAAM. KOA 
940 FORMAT (2/SX.7H0A MAX-,F7.2.f;:;.2HK-. I3/5X.0HPNL MAX-.F7.2.5X.2HK 
CI3/5X.6HC MAX" . F6 . 2 . QX. 2H!< = . I3/5X.9HPHLT MAX-.F7.2.4X.2HK-. 13/ 
CSX. GIISLA MAX= . F7 . 2 . 5X. 2HK ■ .. 13/) 

UR I Td (6. 950 ) DUR. LT 1 1 . LT2. KA V 

950 F0P"'lT(5X. 15MDURATI0N CORR. ".F6.2.2X.6H. K- 13. IX. 3HAND. 14// 
C5X.25MSPECTRA AVG. TECHIIIPUE = .11/ 

CSX. n: II " LIN./5X.0H2 " EXP./SX. 0H3 = NOME///) 

IE (KT0.E0.2) GO TO 952 
UR ITE (6. 951) 


951 FORNAT(5X.31PSEUDOTONE CORRECTIONS II ICl.UDcD, //) 
GO TO 954 


952 UR ITE (6. 953) 

953 FORM 1TC5X.25HII0 P'JEUDOTONZ CORRECTIONS.//) 
254 UR ITE (6. 955 ) EPML 

955 F0miT(5X. 51 IZPIIL ^ . F7 . 2. 4/) 


UR ITE (1.950) 

960 FO.RNATC ADDITIONAL DATA COPY 7 YES-TYPE 1. NO-TYPE 2*) 
READ(1.*)KC0P 
IF ( l-kCOP)959»410 


959 UR ITE (1.951) 

961 FORMAT ( "DATA SPECTRUM PLOT?"/ 

C" 1YPE: 1 "AMBIENT DATA PLOT"/ 

C" 2= INPUT DATA PLOT"/ 

C" 3=N0 PLOTS DESIRED"/) 

READ( l.*)KSPC 
IF (KSPC-2)963.970.992 

970 UR ITE (1.971) 

971 FORMAT ( "ENTER NO. OF SPECTRUI1 TO GE PLOTTED"/ 
C "ENTER 0 IF NO PLOTS DESIRED"/) 

READ( l.*)K 

962 IE (K.EQ. 0)992. 972 
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i 

1 


i 


0152 963 SMAX-0. 

0153 DO 9S5 1-1.24 

0154 SPDO( I) -SPM( I) 

0155 965 SMAX-AMAX1(SPD0(I),SMAX> 

0156 K-0 

0157 GO TO 969 

0150 C 

0159 C SELECTED SPECTRUM PLOT 

0160 C 

0161 972 SMAX-0 

0162 DO 973 1-1.24 

0163 5PD0(I)“DAV(l.K) 

0164 973 SMAX-AMAX1 (SPDO( 1) .SMAX) 

0165 969 IF<AM0D(3MAX. 10.) .GE.5. )974.9?5 

0166 974 SmX-lO.*<AINT(SMAX/'10.) + l.) 

0167 GO TO 976 

0153 975 SMAX- |0.x A IMT (SMAX/ 10. )+5. 

0169 976 SL0-SMAX-30. 

0170 l 1 *? ITE (6. 977)SLQ. K» SMAX 

0171 977 FORMAT <20^6X» 1!I>.X.F3.0.X.2HD0. 15X.2HK-. I3.32X. IH<,X.F3.0.X.2HDB) 

0172 UR1TEC6.978) 

0173 970 F0RMAT(/X.4MDAND.66X.3HSPL/) 

0174 DO 939 1-1.24 

0175 IMS - IF 1X(2.>k(SMAX-SPD0< I ) ) ) 

0176 IF(59-INS)930.9A2 

0177 900 INS-59 

0170 902 DO 904 M= 1.60- INS 

0179 934 ID (M) -S3 

0100 DO 9S5 M»6t- IMS. 62 

0131 906 1D(MJ -32 

0132 UR I TE (6.933) ICXDC I) . ID.SPDO( I) 

0103 903 FORMATC I4.2X.62R1.F7.2) 

0104 909 CONTINUE 

U1G5 1<RITE(6.990> 

0185 990 Foi?NAT(I0s> 

0187 GO TO 970 

0103 C 

01C9 C FINISH PROGRAM OR LOOP BOCK FOR MORE DATA 

0190 C 

0191 992 CALL EPNL3 

0192 C 

0193 C 

0194 C IMPORTANT PARAMETERS DEFINED 

0195 C 

0196 C DC: CAL GAIN 

0197 C BD: DATA GAIN 

0193 C BDA: AMU l ENT NOISE GAIN 

0199 C CMAX: MAX. TONE CORRECTION 

0200 C DAV: AVERAGE DATA SPECTRUM 

0201 C DUR: DURATION CORRECTION 

0202 C OA: OVERALL DATA SPECTRUM 


2L 
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0^33 C 000: OVERflLL *fp tJT. 

02C4 C Oft^Ms OVERP.LL *A" UT. MAX. 

0205 C Or.M: OVERALL MAX. 

O’PS C r Ci PINK NOISE CORRECTION 

CW C PCI.: PERCEIVED NOISE LEVEL 

0203 C Pi LT: TUNE CORRECTED PERCEIVED NOISE LEVEL 

3 c I'M. Til: TONS CORRECTED PERCEIVED NOISE LEVEL 

Ip IQ C PC: MICROPHONE CORRECTION 

0211 C SI'iT : AVERAGE PINK NOISE SPECTRUM 

0212 C SPOIl: rtf 12 1 ENT NOISE SPECTRUM 

0213 C SPD: DATA SPECTRUM 

0214 C STDC : CORRECTED DATA SPECTRUM 

U215 C SPP: PINK NOISE SPECTRUM 

0216 C 

021? C 

0213 399 END 


x-K NO ERRORS* 
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1 TAPE MOTION CONTROLS 

•2 SERIAL TIME CODE AND TEST CONDITION INFORMATION (GAlV, ETC.) 
f 3 ANALOG NOISE DATA 

[\ PARALLEL TIME CODE AND TEST CONDITION INFORM. ION 

5 DIGITAL (PARALLEL) 1/3 OCTAVE FILTER DATA 

6 CONTROL LINES (DATA READY, ADVANCE CHANNEL, ETC.) 

7 a. PARALLEL TIME CODE AND TEST CONDITION INFORMATION 
b. "STAPT TIME" SIGNAL 

6 CONTROL WORDS FOR: 

a. READ TIME INFORMATION 

b. READ TEST CONDITION INFORMATION 

c. MOTION COMMANDS FOR ANALOG TAPE RECORDER 

9 RECORDING OF PROCESSED DATA 

(JL<£ COMMUNICATION LINK FOR OPERATOR 
11 DISC OPERATING SYSTEM COMMUNICATION 


Figure 1.- System block diagram. 



Figure 2.- Minicomputer and standard peripherals (CRT terminal, 
disc, magnetic tape, printer/plotter) 
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Figure h.- One-third octave analyzer. 
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figure 6 .- 


Time 


code 


reader/tape 


search unit. 


START TIME AT K*l- 12:13:59: 4 


RWRODUCIBII.rpy OK thp 
PAOB IS KO(™ 


K 

1 X 

2 XXX 

3 XXX 

4 XXXXXXX 

5 XXX 

6 XXXXX 

7 XXXXX 

8 XXXXXXXXXX 

9 xxxxxxxxxxxxxxxxxxx 

1 0 xxxxxxxxxxxxxxxxxx 

1 1 xxxxxxxxxxxx 

12 xxxxxxxxxxxxxxxxxxxxxx 

1 3 X>900<XXX><XXXXXXXXXXXXX 

14 xxxxxxxxxxxxxxxxxxxxxxxx 

15 xxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

16 xxxxxxxxxxxxxxxxxxxxxxx 

17 xxxxxxxxxxxxxxxxxxxxxx 

18 xxxxxxxxxxxxxxxxxxxxxxxxx 

19 xxxxxxxxxxxxxxxxxxxxxxx 

20 xxxxxxxxxxxxxxxxxxx 

2 1 xxxxxxxxxxxxx 

22 XXXXXXXXXXX>CXX 

23 XXXXXXXXXXXXXX 

24 XXXXXXXXXX.XXXXX 

25 XXXX 

26 XX 

27 XXXXX 
23 XXXXX 

29 X 

30 X 


PNLT 

PNL 

OASPL 

SLA 

95.40 

94.77 

86.30 

82.41 

100.04 

100.04 

90.84 

87.90 

100 . 12 

99.62 

89.99 

87.46 

101.95 

100.53 

91.40 

88.06 

99.99 

98.82 

89.99 

86.21 

100.73 

100.23 

91.39 

88.04 

100.79 

98.95 

91.03 

86.82 

103.60 

101.14 

91.01 

87.95 

108. 17 

105.34 

96.09 

92.49 

107.34 

103.96 

94.40 

90.07 

104.40 

102.79 

94.25 

89. 16 

109.40 

107.85 

96.94 

93.87 

109.55 

107.32 

96.25 

92.20 

1 10.61 

107.99 

97.30 

93.21 

113.72 

110.90 

98.96 

95.93 

110.11 

107.54 

96.65 

92.97 

109.57 

106.96 

97.60 

92.53 

111.16 

108.45 

98.95 

93.91 

110.47 

108.24 

99.44 

94.52 

107.98 

106.59 

99.60 

93.37 

105.02 

104.41 

98.48 

91.83 

105.36 

104.24 

98.73 

91.69 

105.47 

104.43 

99.31 

90.81 

105.90 

104.02 

98.00 

91. 18 

100.23 

98.65 

94.78 

84.81 

99.42 

90.05 

93.51 

83.39 

101 . 16 

99.96 

96.45 

S5.46 

103.82 

99.78 

97.03 

85.51 

96.77 

95.57 

94.62 

80.49 

94.49 

93.24 

92.63 

78 . 03 


0A MAX- 99.60 K- 20 
PNL MAX- 110.90 K» 15 
C MAX- .00 K- 10 
PNLT MAX- 113.72 K- 15 
SLA MAX- 115.71 l<« 13 


DURATION CORR.n -5.53 . K- 9 AND 24 


SPECTRA AVG. TECHHIOUE - 3 

1 - LIN. 

2 - EXP. 

3 - NONE 


NO PSEUDOTONE CORRECTIONS 
EPNL- 100. 18 


Figure 7.- PNLT Plot 



1 Report No 2 Government Act avion No 

NASA TM 7871 3 


4 t itte jno Subtitle 

A I’OMIHV l<- •OiTI’RiH.IiKI' IN. "I'IUIMi. STATION ;".YKTKM HW Til 1 HU- 
iVTAVl' AiiAl.Y^U’ 

1 Authorisl 

jett.’u . : m ‘in and .lame; - . II. Monte ith 


9 Peitufinmg Oigenization Name end Address 

Lai i ,* i ey beaearch Center 
llumpt ii, 1 I Sbb 


I? Sponsoring Agency Name and Address 

lati.iKJ Aeronautics an.i Spare Administration 
.•.ushington, DC 


1 Recipient s Catalog No 


S Report Date 

Apr! I '1 , 1978 


6 Performing Organization Code 

icy 5 


8 Performing Organization Report No 


10 Work Unit No 


1 1 Contract or Grant No 


13 Type of Report and Period Covered 

Technical Memcrandun 

14 Sponsor mg Agency Code 

505-03-1 3-1 ^ 



16 Abstract 

tation tei Is iescribed which employs a minicomputer, a Dne- third 
ictave band analyzer, and a time code/tape search unit for the automatic contr 
and analysis third-octave data. With this system the Information necessary 
for iata adjustment is formatted in such a way as to eliminate much operator 
Interface, thereby substantially reducing the probability for eriar. A 
iescripti n of a program for the calculation of effective perceived noise level 
(KFiJL) fr m aircraft noise lata is included as an example of how this system 
can be use i. 


1/ Key Word* (Suggested by Authorisl) 

Minicomputer 

Time Code/Tape Learch Unit 
KFNL Calculation 


18 Distribution Statement 


Unclassified - Unlimited 


19 Security Claud (of this reportl 

Uric 1 as si fied 


20 Security Claisd (of this page) 

21 No of Pages 

22 Price" 

llncl ass i fied 

32 

$4.50 


’ F° ( sale by the National Technical Information Service Springfield Virginia 22161 


















